APPENDIX A: Source Code for verification of operation of Eye-D 

^include "main.h" 

^define XDIM S12L 
Idefine XDIMR 512 
#define YDIM 480L 
^define BITS 8 
#define RMSVAL 5.0 
#define NUM_NOISY 16 
#define NUM_DEMOS 3 
#define GRAD THRESHOLD 10 

struct char_buf { 

char filename[80]; 
FILE *fp; 
fpost fpos; 
char buf[XDIMR]; 

}; 

struct uchar_buf { 

char filename[80]; 
FILE *fp; 
fpos_t fjpos; 

unsigned char buf[XDIMR]; 

}; 

struct int_buf { 

char filename[80]; 
FILE *fp; 
fpost fpos; 
int buffXDIMR]; 

}; 

struct cortexjs { 

char filename[80]; 
FILE *fp; 
fpos_t fpos; 

unsigned char buf[XDIMR]; 

}; 

struct ucharbuf test_image; 
struct char_buf snow_composite; 
struct ucharjwf distributed_image; 
struct ucharjbuf tempimage; 
struct int_buf tempwordbuffer; 
struct intbuf temp_wordbuffer2; 
struct uchar_buf snow_images; 
struct cortex_s cortex; 

int demo=0; /* which demo is being performed, see notes */ 
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EXPRESS MAIL TB491255248US 



int our_code; /* id value embedded onto image */ 
int found_code=0; /* holder for found code*/ 



int waitvbb(void){ 

while( (_mp(PORT_BASE)&8) ); 
while( !Cinp(PORT_BASE)&8) ); 
return(l); 

} 

int grabb(void){ 

waitvbbO; 

_outp(PORT_BASE+ 1,0); 

_outp(PORT_BASE,8); 

waitvbbO; 

waitvbbO; 

j)utp(PORTJBASE,0xl0); 
return(l); 

} 

int livee(void){ 

j)utp(PORT_BASE,0x00); 
return(l); 

} 



int live_video(void){ 
HveeO; 
return(l); 

} 

int freeze_frame(void){ 
grabbO; 
return(l); 

} 

int grab_frame(struct ucharjmf *image){ 
long i; 

grabbO; 

fsetpos(image- > fp, &image- > fpos ); 
fsetpos(cortex.fp, &cortex.fpos ); 
for(i=0;i<YDIM;i++){ 

fread(cortex.buf,sizeof(unsigned char),XDIMR,cortex.fp); 

fwrite(cortex.buf,sizeof(unsigned char),XDIMR,image-> fp); 

HveeO; 
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return(l); 

} 

int wait_vertical_blanks(int number) { 
long i; 

for(i =0;i < number;i+ +)waitvM>0; 
retura(l); 

} 



int clear_char_image(struct char_buf *charbuffer){ 
long i j; 
char *pchar; 
fpos_t tmp_fpos; 

fsetpos(charbuffer- > fp, &charbuffer- > fpos ); 
for(i=0;i<YDIM;i++){ 

fgetpos(charbuffer->^>, &tmp_fpos ); 

pchar = charbuffer->buf; 

fread(charbuffer- > buf,sizeof(char),XDIMR,charbuffer- > fp); 
for(j =0;j < XDIM;j + +) *(pchar+ +) = 0; , 
fsetpos(charbuffer->fp, &tmp_fpos ); 
fwrite(charbuffer-> buf,sizeof(char),XDIMR,charbuffer- > fp); 

} 

return(l); 

} 

int display_uchar(struct uchar_buf *image,int stretch){ 
unsigned char *pimage; 
unsigned char highest = 0; 
unsigned char lowest = 255; 
long i j; 

double dtemp,scale,dlowest; 
fpos_t tmp_fpos; 

if(stretch){ 

fsetpos(image->fp, &image->fpos ); 

fread(image- > buf,sizeof(unsigned char),XDIMR,image- > fp); 

fread(image-> buf,sizeof(unsigned char),XDIMR, image- > fp); 

for(i=2;i<(YDIM-2);i++){ 

fread(image- > buf,sizeof(unsigned char),XDIMR,image- > fp); 
pimage = &image->buf[3]; 
fora=3;j<(XDIM-3);j + +){ 

if( *pimage > highest )highest = *pimage; 
if( *pimage < lowest )lowest = *pimage; 
pimage++; 

} 

} 
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if(highest = = lowest ){ 

printf("something wrong in contrast stretch, zero contrast"); 
exit(l); 

} 

scale = 255.0 / ( (double)highest - (double)lowest ); 
dlowest = (double)lowest; 
fsetpos(image- > fp, &image- > fpos ); 
for(i=0;i<YDIM;i++){ 

fgetpos(image->fp, &tmpJpos ); 

fread(image- > buf,sizeof(unsigned char),XDIMR,image- > fp); 

pimage = image- >buf; 

forO=0;j<XDIM;j++){ 

dtemp = ((double) *pimage - dlowest)*scale; 

if(dtemp < 0.0)*(pimage++) = 0; 

else if(dtemp > 255.0)*(pimage+ +) = 255; 

else *(pimage-f +) = (unsigned char)dtemp; 

} 

feetpos(image->fp, &tmp_fjpos ); 

fwrite(image- > buf,sizeof(unsigned char),XDIMR,image- > fp); 



fsetpos(image->fp, &image->fpos ); 
fsetpos(cortex.fp, &cortex.fpos ); 
for(i=0;i<YDIM;i++){ 

fread(image- > buf,sizeof(unsigned char),XDIMR,image- > fp); 

fwrite(image- > buf,sizeof(unsigned char),XDIMR,cortex.fp); 

} 

retura(l); 



int clear_int_image(struct int_buf *wordbuffer){ 
long i,j; 
int *pword; 
fpos_ttmp_fpos; 

fsetpos(wordbuffer->fp, &wordbuffer->fpos ); 
for(i=0;i<YDIM;i++){ 

fgetpos(wordbuffer->fp, &tmp Jpos ); 

pword = wordbuffer->buf; 

fread(wordbuffer- > buf,sizeof(int),XDIMR,wordbuffer- > fp); 
for(j=0;j<XDIM;j++) *(pword++) = 0; 
fsetpos(wordbuffer->fp, &tmp_fpos ); 
fwrite(wordbuffer- > buf,sizeof(int),XDIMR,wordbuffer- > fp); 
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retum(l); 

} 



double find mean_int(struct int_buf *wordbufferj{ 
long i,j; 
int *pword; 
double mean =0.0; 

fcetpos(wordbuffer->ip, &wordbuffer->fpos ); 
for(i=0;i<YDIM;i++){ 

pword = wordbuffer->buf; 

fread(wordbuffer- > buf ,sizeof(int),XDIMR,wordbuffer- > fp); 
for(j=0;j<XDIM;j++)mean += (double) *(pword++); 

} 

mean /= ((double)XDIM * (double)YDIM); 
return(mean); 

} 

int add_uchar_to_int(struct ucharbuf *image,struct int_buf *word){ 
unsigned char *pimage; 
int "pword; 
long i,j; 

fpos_t tmpfpos; 

fsetpos(image->fp, &image- > fpos ); 
fsetpos(word- > fp, &word- > fpos ); 
for(i=0;i<YDIM;i++){ 

pword = word- > buf; 

fgetpos(word->fp, &tmp_fpos ); 

fread(word- > buf,sizeof(int),XDIMR, word- > fp); 

pimage = image- > buf; 

fread(image- > buf,sizeof(unsigned char),XDIMR,image- > fp); 
for(j=0;j<XDIM;j + +) *(pword++) += (int)*(pimage++); 
fsetpos(word->fp, &tmp fpos ); 
fwrite(word-> buf, sizeof(int),XDIMR, word- > fp); 

} 

return(l); 

} 



int add_char_to_uchar_creating_uchar(struct charbuf *cimage, 
struct uchar_buf *image, 
struct ucharbuf *out_image){ 
unsigned char *p image, *pout_image; 
char *pcimage; 
int temp; 
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long i j; 

fcetpos(image->fp, &image->fpos ); 
fsetpos(out_image->fp, &out_image->fjpos ); 
fsetpos(cimage->fp, &cimage->fpos ); 
for(i=0;i<YDIM;i++){ 

pcimage = cimage->buf; 

fread(cimage- > buf,sizeof(char),XDIMR,cimage- > fp); 
p image = image- >buf; 

fread(image- > buf,sizeof(unsigned char), XDIMR, image- > fp); 
pout_image = out_image->buf; 
forO=0;j<XDIM;j++){ 

temp = (int) *(pimage++) + (int) *(pcimage++); 

if(temp<0)temp = 0; 

else if(temp > 255)temp = 255; 

*(pout_image++) = (unsigned char)temp; 

} 

fwrite(out_image- > buf,sizeof(unsigned char),XDIMR,out_image- > fp); 

} 

return(l); 

} 



int copy_int_to_int(struct int_buf *word2,struct int_buf *word){ 
long i; 

fsetpos(word2->fp, &word2->fpos ); 
fsetpos(word- > fp, &word- > fpos ); 

for(i=0;i<YDIM;i++){ 

fread(word- > buf ,sizeof(int),XDIMR,word- > fp); 
fwrite(word- > buf,sizeof(int),XDIMR,word2- > fp); 

} 

retura(l); 



void get_snow_images(void){ 

unsigned char *psnow,*ptemp; 
int numberjsnowinputs; 

int temp,*pword,*pword2,bit; 
long i, j; 

double rms,dtemp; 

live_videoO; /* device specific */ 

printf("\n\nPlease point camera at a medium lit blank wall. "); 
printf("\nDefocus the lens a bit as well H ); 
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printf("\nlf possible, place the camera into its highest gain, and "); 
printf("\nput the gamma to 1.0."); 
printf(" Ensure that the video is not saturated "); 
printf("\nPress any key when ready... "); 

while( IkbhitO ); 

printf("\nNow finding difference frame rms value... "); 

/* subtract one image from another, find the rms difference */ 

live_video(); 
wait_vertical_blanks(2); 
grab_frame(&temp_image) ; 
live_video(); 
wait_vertical_blanks(2); 

grab_frame(&distributed_image); /* use first image as buffer */ 
rms - 0.0; 

fsetpos(temp_image.fp, &temp_image.fpos ); 
fsetpos(distributed_image.fp, <Mistributed_image.fpos ); 
for(i=0;i<YDIM;i++){ 

ptemp = tempimage.buf; 

fread(temp_image.buf,sizeof(unsigned char),XDIMR,temp_image.fp); 
psnow = distributed_image.buf; 
fread(distributed_image.buf,sizeof(unsigned 
char),XDIMR,distributed_image.fp); 

forG=0;j<XDIM;j++){ 

temp = (int) *(psnow++) - (int) *(ptemp++); 

dtemp = (double)temp; 

dtemp *= dtemp; 

rms + = dtemp; 

} 

} 

rms /= ( (double)XDIM * (double)YDIM ); 
rms = sqrt(rms); 

printf("\n\nAn rms frame difference noise value of %lf was found.", rms); 
printf("\nWe want at least % If for good measure" ,RMS_VAL); 
/* we want rms to be at least RMS VAL DN, so ... */ 
if(rms > RMS_VAL) number_snow_inputs - 1; 
else { 

dtemp = RMSJVAL / rms; 
dtemp *= dtemp; 

numbersnowjnputs = 1 + (int)dtemp; 

} m 

printf("\n%d images will achieve this noise level", number_snow_inputs); 
/* now create each snowy image */ 
printf("\nStarting to create snow pictures... \n H ); 
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fsetpos(snow_images.fp, &snow_images.fpos ); /* set on first image*/ 
for(bit = 0; bit < BITS; bit+ +){ 

clear_inMmage(&temp_wordbuffer); 
for(i=0;i < number_snow_inputs;i+ +){ 

live_videoO; 

wait_vertical_blanks(2); 

grab_frame(&temp_image); 

add_uchar__to_int(&temp_image,&temp_wordbuffer); 

} 

clear_int_image(&temp_wordbuffer2); 
for(i=0;i < number_snow_inputs;i+ +){ 

live_videoO; 

wait_vertical_blanks(2); 

grab_frame(&temp_image); 

add_uchar_to_int(toemp_image,&temp_wordbuffer2); 

} 

/* now load snow_images[bit] with the difference frame biased by 
128 in an unsigned char form just to keep things clean */ 
/* display it on cortex also */ 

fsetpos(temp_wordbuffer2.fp, &temp_wordbuffer2.fpos ); 
fsetpos(temp_wordbuffer.fp, &temp_wordbuffer.fpos ); 
fsetpos(temp_image.fp, &temp_image.fjpos ); 
for(i=0;i<YDIM;i++){ 

pword = temp_wordbuffer.buf; 

fread(temp_wordbuffer.buf,sizeof(int),XDIMR,temp_wordbuffer.§)); 
pword2 = temp_wordbuffer2.buf; 

fread(temp_wordbuffer2.buf,sizeof(int),XDIMR,temp_wordbuffer2.^)); 

psnow = snow_images.buf; 

ptemp = tempimage.buf; 

fora=0;j<XDIM;j + +){ 

*(psnow++) = *(ptemp++) = (unsigned char) 
(*(pword++) - *(pword2++) + 128); 

} 

fwrite(snow_images.buf,sizeof (unsigned 
char),XDIMR,snowJmages.fp); 

fwrite(temp_image.buf,sizeof(unsigned 
char),XDIMR,temp_image.lp); 

} 

freeze_frameO; 

display_uchar(&temp_image,0); /*1 signifies to stretch the contrast*/ 

printf( H \rDone snowy %d M ,bit); 

wait_vertical_blanks(30); 
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return; 

} 



void loop_visual(void){ 

unsigned char *psnow; 

char *pcomp; 

long i,j, count = 0; 

int ok=0,temp,bit,add_it; 

double scale = 1.0 / RMS_VAL; 

double dtemp,tmpscale; 

fpost tmp_fpos; 

/* initial nns of each snowy image should be around 5 to 10 DN. 
let's assume it is 5, and assume further that our acceptable noise level of 
the full snowy composite should be approximately 1 DN, thus we need to 
scale them down by approximately 5*BITS as a first guess, then do the 
visual loop to zoom in on final acceptable value */ 



printf("\n\n Now calculating initial guess at amplitude... \n w ); 
while( !ok ){ 

/* calculate snow_composite */ 

/* clear composite */ 

clear_char_image(&snow_composite); 

fsetpos(snow_images.fp, &snow_images.fpos ); /* set on first image*/ 
for(bit=0;bit<BITS;bit+ +){ 

j = 128 >> bit; 

if( ourcode & j)add_it= 1; 

else add_it=0; 

fsetpos(snow_composite.lp, &snow_composite.fpos ); 
for(i=0;i<YDIM;i++){ 

psnow = snow_images.buf; 

fread(snow_images.buf,sizeof(unsigned 

char),XDIMR,snow_images.fp); 

fgetpos(snow_composite.fp, &tmp_fpos ); 

fread(snow_composite.buf,sizeof(char),XDIMR,snow_composite.fp); 

pcomp = snow_composite.buf; 
fora=0;j<XDIM;j-h+) { 

dtemp = ((double)*(psnow++) -128.0) * scale; 
if(dtemp<0.0){ 

temp = -(int) fabs( -dtemp +0.5); 
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) 

else { 

temp = (int) fabs( dtemp +0.5); 

} 

if(temp > 127) { 

temp = 127; 

} 

else if(temp < -128) { 
temp = -128; 

} 

if(add_it){ 

*(pcomp++) += (char)temp; 

} 

else { 

*(pcomp+ +) -= (char)temp; 

} 

fsetpos(snow_composite.fp, &tmp_fpos ); 
fwrite(snowjx>mposite.buf,sizeof(char),XDIN^ 

} 

printf("\rDone snowy %d ",bit); 

} 

/* add snow composite to test image to form dist image */ 
add_char_to_uchar_creating_uchar( 

&snow_composite, 

&test_image, 

&distributed_image); 

/* display both and cue for putting scale down, up or ok */ 
i= count = 0; 

printf("\n Depress any key to toggle, enter to move on...\n H ); 
printf("\r Distributed Image... "); 
display_uchar(&distributed_image,0); 
while(getchO != '\r* ){^ 

if( (count* +) % 2){ 

printf("\r Distributed Image. . . "); 

display_uchar(&distributed_image,0); 

} 

else { 

printf("\r Original Image... "); 
display_uchar(&test_image,0); 

} 

} 

printfOnScale = %lf ".scale); 
printf("\nEnter new scale, or > le6 for ok... "); 
scanf(" %ir,&tmpscale); 
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} 



if(tmpscale > le6)ok=l; 
else scale = tmpscale; 

} 

/* distributed image now is ok; calculate actual snow_images used and 
store in those arrays; */ 

fsetpos(snow_images.fp, &snow_images.fpos ); /* set on first image*/ 
printf("\nNow storing snow images as used... \n"); 
for(bit=0;bit<BITS;bit+ +){ 

for(i=0;i<YDIM;i++){ 

psnow = snow_images.buf; 
fgetpos(snow_images.fp, &tmp_fpos ); 
fread(snow_images.buf,sizeof(unsigned 
char),XDIMR,snow_images.fp); 

forO =0;j < XDIM;j + +) { 

dtemp = ((double)*psnow -128.0) * scale; 
if(dtemp<0.0){ 

temp = -(int) fabs( -dtemp +0.5); 

} 

else { 

temp = (int) fabs( dtemp +0.5); 

} 

*(psnow+ +) = (unsigned char)(temp + 128); 

} 

fsetpos(snow_images.fp, &tmp_fpos ); 
fwrite(snow_images.buf,sizeof(unsigned 
char),XDIMR,snow_images.fp); 

} 

printf("\rDone snowy %d ",bit); 

} 

return; 

} 



double findjjrad(struct intjbuf *image,int load_buffer2){ 
int bufl[XDIMR],buf2[XDIMR],buGpCDIMR]; 
int *pbufl,*pbuf2,*pbuf3,*p2; 
double total =0.0,dtemp; 
long i, j; 

fpos_t tmp_pos; 

fsetpos(image- > fp, &image- > fjpos ); 
fgetpos(image->fj3, &tmpjx>s ); 
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fsetpos(temp_wordbuffer2.fp, &temp_wordbuffer2.fpos ); 

for(i= l;i < (YDIM-l);i+ +){ 

fsetpos(image->fp, &tmp_pos ); 
fread(buf 1 ,sizeof(int),XDIMR,image- > fp); 
fgetpos(image->fp, &tmp_pos ); 
fread(buf2,sizeof(int),XDIMR,image- > fp); 
fread(buf3,sizeof(int),XDIMR,image- > fp); 
pbufl=bufl; 
pbuf2=buf2; 
pbu£3=buf3; 

p2 = temp_wordbuffer2.buf; 

if(load_buffer2){ 

for(j=l;j<(XDIM-l);j++){ 

dtemp = (double)*(pbufl + +); 
dtemp += (double)*(pbufl + +); 
dtemp += (double)*(pbufl-); 
dtemp += (double)*(pbuf2++); 
dtemp -= (8.0 * (double) *(pbuf2++)); 
dtemp += (double)*(pbuf2-); 
dtemp += (double)*(pbuf3++); 
dtemp += (double)*(pbuf3++); 
dtemp += (double)*(pbuf3-); 
*p2 = (int)dtemp; 
if( *p2 > GRAD_THRESHOLD ){ 

*(p2++)-= GRAD_THRESHOLD; 

} 

else if( *p2 < -GRAD_THRESHOLD ){ 

*(p2 + +) += GRAD_THRESHOLD; 

} 

else { 

*(p2++) = 0; 

} 

} 

fwrite(temp_wordbuffer2.buf,sizeof0nt) > XDIMR,temp_wordbuffer2.fp); 

} 

else { 

fread(temp_wordbuffer2.buf,sizeof(int),XDIMR,temp_wordbuffer2.fp); 
forO=ly<(XDIM-l);j++){ 

dtemp = (double)*(pbufl + +); 

dtemp += (doubIe)*(pbufl + +); 

dtemp + = (double)*(pbufl-); 

dtemp += (double)*(pbuf2+ +); 

dtemp -= (8.0 * (double) *(pbu£2++)); 

dtemp += (double)*(pbuf2-); 

dtemp += (double)*(pbuf3++); 

dtemp += (double)*(pbuf3++); 
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} 
} 

retura(total); 



dtemp + = (doubIe)*(pbuf3~); 

dtemp -= (double) *(p2++); 

dtemp *= dtemp; 
total + = dtemp; 



void searchl (struct ucharjmf *suspect){ 
unsigned char *psuspect,*psnow; 
int bit,*pword,temp; 
long i j; 

double add_metric,subtract_metric; 
fposj tmp_fpos; 



/* this algorithm is conceptually the simplest. The idea is to step 
through each bit at a time and merely see if adding or subtracting the 
individual snowy picture minimizes some 'contrast' metric. 
This should be the most crude and inefficient, no where to go but 
better */ 



fsetpos(snow_images.fp, &snow_images.fpos ); 
temp =256; 

clear_int_image(&temp_wordbuffer); 
add_uchar_to_int(suspect,&temp_wordbuffer); 

find _grad(&temp_wordbuffer,l); /* 1 means load temp_wordbuffer2 */ 
for(bit=0;bit<BITS;bit+ +){ 
/* add first */ 

fgetpos(snow_images.fp, &tmpJpos ); 
fsetpos(suspect- > lp, Asuspect- > lpos ); 
fsetpos(temp_wordbuffer.f£, &temp_wordbuffer.fpos ); 

for(i=0;i<YDIM;i++){ 

pword = tempwordbuffer.buf; 
psuspect = suspect- >buf; 
psnow = snowimages.buf; 

fread(suspect->buf,sizeof(unsigned char),XDIMR,suspect->fp); 
fread(snow_images,buf,sizeof(unsigned 
char),XDIMR,snow_images.fp); 

forG=0;j<XDIM;j + +){ 

*(pword+ +)=(int)*(psuspect+ +)+(int)*(psnow+ +)-128; 

} 
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fwrite(temp_wordbuffer.buf,sizeof(int),XDIMR,temp_wordbuffer.fp); 

} 

add_metric = find_grad(&temp_wordbuffer,0); 
/* then subtract */ 

fsetpos(snow_images.fp, &tmp_fpos ); 
fsetpos(suspect->^>, &suspect->fpos ); 
fsetpos(temp_wordbuffer.fjp, &temp_wordbuffer.fpos ); 
for(i=0;i<YDIM;i++){ 

pword = tempwordbuffer.buf; 

psuspect = suspect- >buf; 

psnow = snowimages.buf; 

fread(suspect- > buf,sizeof(unsigned char), XDIMR, suspect- > fp); 
fread(snow_images.buf,sizeof(unsigned 
char),XDIMR,snow_images.fjp); 

forQ=0;j<XDIM;j++){ 

*(pword+ +)=(int)*ft)suspect+ +)-(int)*(psnow+ 4-)+ 128; 

} 

iwrite(temp_wordbuffer.buf,s^ 

} 

subtractmetric = find_grad(&temp_wordbuffer,0); 

printf( w \nbit place %d: add=%le , 
sub = % le" ,bii,add_metric,subtract_metric); 
temp/ =2; 

if(add_metric < subtract_metric){ 
printff bit value = 0"); 

} 

else { 

printff bit value = 1 H ); 
foundcode + = temp; 

} 

} 

printf("\n\nYour magic number was %d",found_code); 
return; 

} 



void search_2(unsigned char *suspect){ 
if(suspect); 
return; 

} 
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void loop_simulation(void){ 

unsigned char *ptemp,*pdist; 
int *pword,int_mean,ok=0,temp; 
long ij; 

double mean,scale; 

/* grab a noisy image into one of the temp buffers */ 
printf("\ngrabbing noisy frame... \n"); 
clear_int_image(&temp_wordbuffer); 
for(i=0;i<NUM_NOISY;i+ +){ 

live_yideo(); 

wait_vertical_blanks(2); 

grab_firame(&temp_image); 

add_uchar_to_int(&temp_image,&temp_wordbuffer); 

j = (Iong)NUM_NOISY; 

printf( H \r%ld of %ld ",i+l j); 

} 

/* find mean value of temp_wordbuffer */ 
mean = find_mean_int (&temp_wordbuff er) ; 
intmean = (int)mean; 

/* now we will add scaled version of this 'corruption' to our distributed 
image */ 
scale = 1.0; 
while( !ok ){ 

/* add noise to dist image storing in temp_image */ 
fsetpos(distributedJmage.fp, &distributed_image.fpos ); 
fsetpos(temp_wordbuffer.fp, &temp_wordbuffer.f]X)s ); 
fsetpos(temp_image.fp, &temp_image.fpos ); 
for(i=0;i<YDIM;i++){ 

pdist = distributed_image.buf; 
pword = temp_wordbuffer.buf; 
ptemp = tempimage.buf; 
fread(distributed_image.buf,sizeof(unsigned 
char),XDIMR,distributed_image.fp); 

fread(temp_wordbuffer.buf,sizeof^^ 
forO=0;j<XDIM;j++){ 

temp = (int) *(pdist++) + *(pword++) - int_mean; 

if(temp<0)temp = 0; 

else if(temp > 255)temp = 255; 

*(ptemp++) = (unsigned char)temp; 

} 

fwrite(temp_image.buf,sizeof(unsigned 
char),XDIMR,temp_image.fp); 

} 



37 




1 ■:>;■*- 'jo^i>{iooU{(i». 



/* display the dist image and the corrupted image */ 
display_uchar(&tempjmage,0); 

/* apply new •corrupted* image to search algorithm 1 for id value */ 
search_l (&temp_image); 

/* apply new 'corrupted' image to search algorithm 2 for id value */ 
/* 

search_2(temp_image); 
*/ 

/* prompt for upping noise content or ok */ 
ok = 1; 

} 

return; 



int initialize_everything(void){ 
long i,j; 

unsigned char *pucbuf; 
char *pcbuf; 
int *pibuf; 

/* initialize cortex */ 
strcpy(cortex.filename,"f:image"); 
if((cortex.fp = fopen(cortex.filename, "rb")) = = NULL){ 
systemfv f g"); 

} 

else fclose(cortex.fp); 

if( (_inp(PORT_BASE) = = OxFF) ){ 

printf("oops "); 

exit(0); 

} 

/* open cortex for read and write */ 

if ((cortex, fp = fopen(cortex. filename, "rb + "))= = NULL){ 

printf( tt No good on open file joe M ); 

exit(0); 

} 

fgetpos(cortex.fp, &cortex.fpos ); 

/* testjmage; original image */ 

st^cpy(tesMmage.filename/e:tstJmg ,, ); 
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if((test_image.fp = fopen(test_image. filename, "wb")) = = NULL){ 
printf(" No good on open file joe "); 
exit(0); 

} 

pucbuf = test_image.buf; 
for(i=0;i<XDIM;i+ +)*(pucbuf+ +)=0; 
for(i=0;i<YDIM;i++){ 

fwrite(test_image.buf,sizeof(unsigned char),XDIMR,test_image.fp); 

} 

fclose(test_image.fp); 

if((test_image.^=fopen(test_image.filename,"rb+ "))= =NULL){ 
printf(" No good on open file joe 
exit(0); 

} 

fgetpos(test_image.fp, &test_image.fpos ); 

/* snow_composite; ultimate image added to original image */ 
strcpy(snowj»mposite.filename,"e:snw_cmp"); 
if((snow_composite.fp=fopen(snow_ro^ 

printfC No good on open file joe H ); 

exit(0); 

} 

pcbuf = snow_composite.buf; 
for(i=0;i<XDIM;i+ +)*(pcbuf+ +)=0; 
for(i=0;i<YDIM;i++){ 

fwrite(snow_composite.buf,sizeof(char),XDINIR,snow_composite.fp); 

} 

fclose(snow_composite.fp); 

if((snow_composite.fp=fopen(snow_composite .filename, "rb+ "))= =NULL){ 
printf(" No good on open file joe "); 
exit(0); 

} 

fgetpos(snow_composite.fp, &snow_composite.Q>os ); 

/* distributedjmage; test_img plus snow_composite */ 
strcpy(distributed_image.filename, H e:dst_img"); 

if((distributed Jmage.fp = fopen(distributed_image. filename, M wb M )) = = NULL){ 
printff No good on open file joe w ); 
exit(0); 

} 

pucbuf = distributed_image.buf; 
for(i=0;i<XDIM;i+ +)*(pucbuf+ +)=0; 
for(i=0;i<YDIM;i++){ 

fwrite(distributed_image.buf,sizeof(unsigned 
char),XDIMR,distributed_image.fp); 

} 

fclose(distributed_image.fp); 

if((distributed_image.fp =fopen(distributed Jmage.filename, "rb 4- H )) = = NULL){ 
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printf(" No good on open file joe "); 
exit(0); 

} 

fgetpos(distributedJmage.fp, &distributed_image.fpos ); 



/* temp_image; buffer if needed */ 

strcpy(tempJmage.filename/e:temp_img H ); 

if((tempJmage.fp=fopen(temp_im^ 

printf(" No good on open file joe "); 
exit(0); 

} 

pucbuf = temp_image.buf; 
for(i=0;i<XDIM;i+ +)*(pucbuf+ +)=0; 
for(i=0;i<YDIM;i++){ 

fwrite(temp_image.buf,sizeof(unsigned char),XDIMR,temp_image.fp); 

} 

fclose(temp_image.fjp); 

if((temp_image.fp = fopen(temp_image.filename, "rb + "))= = NULL){ 
printf( H No good on open file joe "); 
exit(0); 

} 

fgetpos(temp_image.fp, &temp_image.fpos ); 

/* temp_wordbuf¥er; 16 bit image buffer for averaging */ 
strq>y(temp_wordbuffer.filename, ,, e:temp_wrd"); 

if((temp_wordbuffer.fp = fopen(temp_wordbuffer.filename, " wb")) = = NULL){ 
printf(" No good on open file joe "); 
exit(0); 

} 

pibuf = temp_wordbuffer.buf; 

for(i=0;i<XDIM;i++)*(pibuf++)=0; 

for(i=0;i<YDIM;i++){ 

fwrite(tempjvordbufTer.buf,ste^ 
fclose(temp_wordbuffer.fp); 

if((temp_wordbuffer.fjp=fopen(temp_wordbuffer,fUename/rb+ =NULL){ 
printf( M No good on open file joe "); 
exit(0); 

} 

fgetpos(temp_wordbuffer.fp, &temp_wordbuffer.fpos ); 

/* temp_wordbuffer2; /* 16 bit image buffer for averaging */ 
strcpy(temp_wo^dbufTer2.filename, r, e:tmp_wrd2 ,, ); 

if((temp_wordbuffer2.fp = fopen(temp_wordbuffer2 .filename, "wb ")) = = NULL){ 
printff No good on open file joe w ); 
exit(0); 

} 
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pibuf = temp_wordbuffer2.buf; 
for(i=0;i<XDIM;i+ +)*(pibuf+ +)=0; 

for(i==0;i<YDIM;i++){ 

fwrite(temp_wordbuffer2.buf,sizeof(int),XDIMR,te^ 

fclose(temp_wordbuffer2.fp); 

if((temp_wordbuffer2.^=fopen(temp_wordbuffer2.filenam =NULL){ 
printf(" No good on open file joe "); 
exit(0); 

} 

fgetpos(temp_wordbuffer2.fp, toemp_wordbuffer2.fpos ); 

/* snow images; BITS number of constituent snowy pictures */ 
strcpy(snow_images.filename,"snw_imgs"); 

if((snow_images.fp = fopen(snow_images.filename, " wb ")) = = NULL){ 
printf(" No good on open file joe "); 
exit(0); 

} 

pucbuf = snowimages.buf; 

for(i=0;i< XDIM;i+ +)*(pucbuf+ +)=0; 

forO=0;j<BITSy++){ 

for(i=0;i<YDIM;i++){ 

fwrite(snow_images.buf,sizeof(unsigned(±ar),XDIMR,snow_images.Q)); 

} 

fclose(snow_images.fp); 

if((snow_images.fp=fopen(snow_images.filename, H rb + "))= =NULL){ 
printf(" No good on open file joe "); 
exit(0); 

} 

fgetpos(snow_images.lp, &snow images. fpos ); 
return(l); 



int close_everything(void){ 

fclose(test_image.lp); 

fclose(snow_composite.fp); 

fclose(distributed_image.^)); 

fclose(temp_image.fp); 

fclose(temp_wordbuffer.fp); 

fclose(temp_wordbuffer2.fi)); 

fclose(sno w_images . fp) ; 

return(l); 

} 
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main(){ 

int i j; 

printf( H \nInitializing...\n\n"); 

initialize_everythingO; /* device specific and global mallocs */ 
Hve_videoO; 

/* prompt for which of the three demos to perform */ 
while( demo < 1 1 1 demo > NUM_DEMOS){ 

printffWhich demo do you want to run?\n\n"); 

printf(" 1 : Digital Imagery and Very High End Photography Simulation\n n ); 
printf("2: Pre-exposed Print Paper and other Dupping\n"); 
printf("3: Pre-exposed Original Film (i.e. In-Camera)\n"); 
printf("\nEnter number and return: * "); 
scanf("%d\&demo); 

if(demo < 1 1 1 demo > NUM_DEMOS){ 
printf( w \n eh eh n ); 

} 

} 

/* acquire test image */ 

printf("\nPress any key after your test scene is ready... "); 
getchO; 

grab_frame(&test_image); /*grab_frame takes care of device specific stuff*/ 

/* prompt for id number, 0 through 255 */ 

printf("\nEnter any number between 0 and 255.\n"); 

printf( M This will be the unique magic code placed into the image: "); 

scanff %d\&our_code); 

while(our_code< 1 || our_code>256){ 

printf( w Between 0 and 255 please 

scanf( w %d H ,&our_code); 

} 

/* feed back the binary code which will be embedded in the image */ 

printf("\nThe binary sequence "); 

for(i=0;i<BITS;i++){ 

j = 128 > > i; 

if( our_code & j)printf( B 1 "); 

else printfCO"); 

} 

printf(" (%d) will be embedded on the image\n",our_code); 

/* now generate the individual snow images */ 
get_snow_images(); 
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loop_visual(); /* this gives visual feedback on 'tolerable* noise level */ 

printf("\nWe're now to the simulated suspect... \n"); 
loop_simulationO; 

close_everythingO; 

retum(0); 

} 
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